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ABSTRACT 


The time-dependent mixed-layer depth and temperature 
response of a two-layer ocean model to a moving hurricane 
model was investigated. The hurricane model was that of 
Elsberry, Pearson, and Corgnati (1974). The maximum hurri- 
cane wind was bounded at 50 m/sec in order to isolate ocean 
reactions from hurricane variance. The ocean model was 
based on that of Fraim (1973) which was derived from the 
earlier mixed-layer models of Kraus and Turner (1967) and 
Denman (1973). The ocean model was an idealized, quasi- 
Mirec-Gimenslonal model designed to examine the effects 
along the path of movement due to the passage of a moving 
hurricane. Results indicated that slower-moving storms 
engendered the greatest change in mixed-layer depth and 
memperature. Away from the center a combination of entrain- 
ment mixing and convection was responsible for the majority 
of the depth and temperature changes. Near the eye, in the 
Upwelling zone, a combination of entrainment mixing, convec- 
tion, and advection caused temperature change while advec- 


tion alone dominated layer depth change. 
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I. INTRODUCTION 


There exists a large body of evidence which demonstrates 
that there is a complex interaction between a hurricane and 
the ocean. For example it appears that sea-surface tempera- 
mimes im excess Of 26 © are required to produce hurricane- 
force winds. Through these winds, the hurricane imparts 
stress to the ocean, which acts as a sink for the hurricane's 
momentum. Interactions have been modeled in several ways. 
/meeis apparent that for an ocean-interacting hurricane model 
to be valid, the heat-momentum energy-exchange balance must 
be satisfied. 

Elsberry, Pearson, and Corgnati (1974) developed an axi- 
Symmetric, ocean-interacting hurricane model which met the 
energy balance conditions. A description of the air-sea 
interaction within the boundary-layer of a hurricane was the 
output of and the primary purpose for the design of the model. 
meretciis purpose, a stationary, mature hurricane over a con- 
Stant, isothermal ocean was assumed. The capability for in- 
teraction with a time-dependent ocean model was included in 
ime design. 

One of the uses envisioned for the model described above 
was to investigate the cause of low sea-surface temperatures 
.found in the wake of a hurricane. Based on limited observa- 
tions following tropical cyclones, wind mixing, convective 
mixing, and upwelling appear to be the main cooling mechan- 


isms. Jordan (1964) concluded that mixing was dominant. 





Fischer (1958) mentioned upwelling as being important. 
Leipper (1967) found upwelling to be dominant near the eye 
with a transition to convective and/or mechanical mixing 
dominance at greater radii. The classical model by O'Brien 
and Reid (1967) and O'Brien (1967) simulated a stationary 
hurricane with constant forcing over a one-layer ocean. Re- 
sults of the model emphasized the importance of upwelling 
and downwelling. Gilbert (1974) modeled a stationary, time- 
dependent hurricane over a multi-layer ocean. Results of 
PiowGilbert model indicate transitory inertial oscillations 
exert as much influence as the steady-state wind. 

Fraim (1973) developed a time-dependent, two-layer model 
for the ocean only based on the earlier mixed-layer models 
of Kraus and Turner (1967) and Denman (1973). The model was 
driven by the time-dependent, axi-symmetric, stationary 
hurricane model of Elsberry, et. al. (1974). The axi-sym- 
metric assumption permitted representation of the hurricane 
and the ocean by a radial cross-section assumed to be the 
integrated average of all radial cross-sections. Fraim 
(1973) concluded that a combination of entrainment mixing 
and upwelling was the major factor in changing the mixed- 
iayer temperature and depth. 

One of several logical extensions of the model was to 
use it to examine the effects on the ocean induced by a 
moving hurricane. One purpose for moving the storm was to 
consider Geisler's (1970) internal wake theory. The theory 


involves a comparison of the hurricane translational speed 
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(U) with the baroclinic wave speed (C). For a slow moving 
storm (U/C < Y2), Geisler postulated upwelling and resultant 
surface cooling. For a fast moving storm (U/C > 2), Geisler 
postulated an internal wave wake at the interface between 
the mixed layer and thermocline. It should be emphasized 
that while Geisler's model and a similar model by O'Brien 
(1968) were primarily hydrodynamic, the model under study 
here was thermodynamic in nature. 

The present model was originally designed to be fully a 
three-dimensional model. However, due to computer time and 
storage constraints, the model was made to represent a ver- 
tical cross-section of the ocean along which the storm was 
moved. The storm was translated at several typical speeds 
(including U = 0) in order to determine the effects of move- 


ment on changes in the mixed-layer depth and temperature. 


hal 





PiteewrocktPlPON OF [THE MODEL 


A. THE HURRICANE MODEL 

iicmivreteanesmodelswas a two-layer, time-dependent model 
by Elsberry, Pearson, and Corgnati (1974). By assuming the 
hurricane to be mature, an assumption of axi-symmetry could 
be made. A surface inflow layer and an upper outflow layer 
constituted the model. The tangential wind (Vo) profile at 


fne top of the inflow layer was taken to be 


1 
Vor = constant (1) 


at radii outside the radius of maximum wind T.. Inside Tr. 
the winds were assumed to obey 


V ri = constant. (2) 


Equations (1) and (2) were empirically derived from observa- 
tions by Riehl (1963) and supported by observations by Gray 
and Shea (1973). The maximum extent of the storm was de- 
fined to be the radius (rQ) at which flow in the outflow 
Miaereturnea irom Gyclonic to anti-cyclonic, 1.e. Vo = Q. 
The radius of maximum wind was assumed to coincide with the 
eyewall radius, which was taken to be the radius at which 
the average of five adjacent values of equivalent potential 
temperature on) was a maximum. Since the eye-wall has been 
observed to be a region of massive latent heat flux due to 


evaporation, the maximum oe was strongly tied to the sea- 


surface temperature of the underlying water. An expression 


Ike 





for the maximum wind speed was determined by assuming a 


gradient balance between the wind field and the pressure 


field 
y2 
Tapes 
fier 5 


(where f is the coriolis parameter) and using the relation 


Vo = U : (4) 





Matecrating equation (3) from r, to r. gives 





02 5x 
V = S [-2f/r. rr. - r- 
G. ae sO a 
C5) 
2 2 6 Pome 
Vey TRL) r.) + 10 “pr, “e, - Be T° 


From r, to 2r, the wind profiie in (1) was modified to de- 

crease exponentially to 5 meters/second as a transition zone 

Bemrealistic boundary conditions in order to allow downwelling. 
From observations by Riehl (1963), the surface pressure 


at the eyewall was defined by the relation 


6p = 2.5666, (6) 


where the 6 symbols indicate a deviation from p = 1005 milli- 
bars and 6,=3550 K. The radial pressure field was computed by 
mieegrating the gradient wind over the radius, using the wind 
meeeda from (1) and (2). 

The atmospheric boundary-layer model consisted of a sur- 
face layer and a spiral layer above. Specification of the 


surface layer was by Monin-Obukhov similarity theory. Within 
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the surface layer, the surface roughness was a function of 
the stress at the surface. The spiral layer, a transition 
layer from radial to tangential flow, was specified by a 
modification of Ekman-layer theory. An iterative solution 
was used to match the two layers. Further details of the 
iterative solution can be found in Cardone (1969). 

A simple model for the hurricane was required to keep 
computer time and storage requirements to a minimum. Compar- 
atively simple specification of the boundary layer and lack 
Pemopecitication of upper layers enabled compliance with 
those requirements. The axi-symmetric assumption allowed 
fmepresentation of the storm by a radial cross-section. Again, 
computational considerations required such an assumption. A 
#riction velocity for the computation of water velocities 
and kinetic energy gain, and the heat flux were provided to 
the ocean model by the hurricane model. A new radial distri- 
bution of mixed-layer temperature was predicted by the ocean 
fedel, and then provided to the hurricane model. This type 
of interaction can lead to storm intensification. However, 
for the purposes of this study, the hurricane maximum wind 
was bounded at 50 m/sec. This enabled isolation of changes 


in the ocean. 


Bee LHE OCEAN MODEL 


Meee Descriperon of the Model 
The ocean model was a modification of that by Fraim 


(1973), which was based on earlier mixed-layer models by 


Kraus and Turner (1967) and Denman (1972). The ocean was 
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assumed to have an initial mixed-layer depth of 30 m with a 
mixed-layer temperature of 30 C. The below-layer gradient 
was initially 0.1 C/meter to the top of an undisturbed layer 
defined by T = 20 C at Z = 130 meters. These values are 
ical of those found in the Gulf of Mexico and the Caribbean 
Sea in June, July, and August. The model was idealistic in 
that truly isothermal mixed layers and linear below-layer 
Bradients are seldom found. Also, it would not be expected 
that the same thermal structure would exist over the whole 
Seed, Or Even any sizeable fraction of it. 

There were several important modifications. of the. 
Fraim (1973) model. The Fraim model examined a radial cross- 
section of the ocean beneath a stationary hurricane. Equa- 
Brons f£Oor the ocean model were written in R-68 co-ordinates. 
imiteepresent model examined a vertical cross-section of the 
@ecan along which the diametrical cross-section of the hur- 
mredne was allowed to move. The diametrical cross-section 
@emsisted of the radial cross-section used by Fraim and the 
mirror image of that radial cross-section. This was made 
possible by the assumption of axi-symmetry in the hurricane 
Medel, tthe ocean equations were written in cartesian co- 
mmainates on a 200 X 3 grid (see Figure 1) with a 6 kilo- 
Smecer grid whereas Fraim used a 100 X 1 grid with a 3 
maeometer prid spacing. The time increment was 8 minutes 
sompared to the 10 minute increment used by Fraim. A time 
increment of 8 minutes was used because it gave finer con- 


meal over available translational speeds. The movement of 


ts 
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the storm was achieved by moving the pentey of the storm to 
successive grid-points at time increments required to ap- 
Proximate the desired speed. Another modification was the 
addition of a horizontal return flow within the thermocline 
layer extending down to a specified layer of no motion. This 
closed the ocean system except for the radial advection, and 
thus improved the heat balance calculation for the stationary 
storm. Lastly, there was added to the model a convective 
Bajustment for local convective instabilities. 

Several important assumptions were made about the 
manner in which the ocean and hurricane would interact. 
First, since the storm was assumed to be mature at initia- 
tion, the wind stress was imparted to the ocean exponential- 
ly so that after 9 hours the ocean felt 95% of the stress. 
Second, some assumption had to be made as to what fraction 
Meme Stress would be used in inducing turbulence for en- 
trainment mixing. The balance of the stress would be used 
miecne Creation of horizontal currents. With no pertinent 
data available it was decided to let entrainment dominate 
Mmacially, An €xXponential transition dependent on the mixed- 
layer depth shifted the dominance to currents as the mixed 
layer deepened. Third, it was assumed that ocean currents 
began in Ekman balance and adjusted immediately to changes 
mete Stress pattern. This assumption eliminates inertial 
waves which must be resolved using smaller time steps. 
Fourth, in examining only a vertical cross-section it was 


assumed that there was no variation in temperature or ocean 
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velocity across the path of the hurricane. This simplifica- 
tion was made because the net Ekman current should be normal 
to wind flowing inward at a specified angle, rather than 
being directed radially. Because of this, there was a re- 
sultant lack of symmetry across the path. This lack of sym- 
Metry introduced instabilities into the model. Fifth, it 
was assumed that the storm was wind dominated, implying that 
downward heat flux due to solar radiation was negligible. 
iiiemw@ccanmnoadel was forced by kinetic energy input 
due to wind stress, upward sensible heat fluxes, and upward 
Maecent heat tilux due to evaporation. Horizontal currents 
due to wind stress were computed. Resultant upwelling and 
downwelling due to divergent and convergent current regimes 
were then determined. A new mixed-layer depth was calculated 
morough kinetic energy, heat flux, and advection considera- 
means. based on the new mixed-layer depth, a new mixed-layer 
Mmaipctatire was GCdlculated. The below-layer gradient and 
Meaanperature at the top of the thermocline were then deter- 
mined. In this model the mixed-layer temperature and the 
fmeipetarire at the top of the thermocline can be different. 
meethe temperature at the top of the thermocline was greater 
than the mixed-layer temperature a convective overturning 
process was initiated in order to stabilize the thermal 
Seoecture. Following these calculations a radially averaged 
-mixed-layer temperature was passed to the radial cross-section 


of the hurricane. 
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2. Derivation of the Model Equations 


Following Denman (1973), the first law of thermo- 


dynamics can be written 





ie 
a (7) 
dt Pwon PWD 
where Qe = Rue Vee tie wneat Source term, 
@. = oGean water density, 


Y = average extinction coefficient, 
mrespeciticieat at Constant pressure, 
Lee-ecolar lactation ancident on the sea-surface, and 
z = depth in the ocean, positive downward. 
Expanding the total derivative in (7) yields 


Q : : 
ieee ot. ol eget (8) 


PWD at 


where w = oz 1s positive for downward motion. 


Time-averaging the turbulent form of (8) gives 





























Rr pile: vT wT me ny 
_ oT , out , ovl , awl , owl’ _ dT. oaw'T' 
P PEM) ys 02) 2) 6 dtCOCOz (9) 
wp 
Or 
aT = fro TT are 
at Owen OZ 


where the terms u'T' and v'T' were neglected because it was 
assumed that there was no systematic correlation between u' 
ama v' and T'. 

Tug noOwe LO an €Xamination of the boundary condi- 


mons, firstly at the ocean surface, z = 0, a heat loss at 
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the sea surface will induce convective overturning which 
cools the mixed layer. Therefore, the downward turbulent 
heat flux wt must equal the non-radiative heat loss, 


ie . 
Ca as (eine 


where Qn - Qe + Qo » 


and Qe Upvaro Latent heat flux due to evaporation, 

and Qc = upward sensible heat flux. 

At the base of the mixed layer downward turbulent heat flux 
mecaused by wind-induced mixing and convective overturning. 
Simemtlroulent heat flux must equal the heat transfer due to 
entrainment of water into the mixed layer. This balance is 


expressed by 


dh 


GT), = Cae 7 w UCTS > TH) (12) 


where h, = Maxedelayer depth, 
T_ = mixed-layer temperature, 
Th -MEchMemAeuremamucmle FCOpeot the thermocline, and 
= aia the time variation of the mixed-layer due to wind 
mixing and convective overturning. 
For entrainment to occur (dh /dt) PRIMUS there can be 
mon shallowing of the mixed layer unless w is upward. Inte- 
grating (10) over the mixed layer and substituting the bound- 
ary conditions gives 


Ae Qr aac Qa 
eemiemecgeam > 9 © de “TS - Ty) - a 
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Peandineg the total derivatives, neglecting solar radiation, 


and solving for oh /at gives 








oh -h aye duT ann owT 
O _ O Sn Se =a =] 
ot Cea.) dt OX dy OZ 
(14) 
: on oh dh, 


; Pycp (1, -Th) ete oe V oy 
which is a predictive equation for deepening shallowing of 
the mixed layer in terms of the local change in mixed-layer 
temperature. | 

To derive a similar predictive equation for the 
mixed-layer temperature Denman (1973) integrated (10) to an 
meabttrary depth z, and then over the mixed-layer depth 


yielding 


ne 0 Q,h 
so S=-f WT, a - 4 


O Meee 





(15) 


h 
where -S ° Gm). G@zetsetne Conversion Of potential energy 
Demkinetic energy by convection. If the turbulent kinetic 
energy 1s not varying in time, we may write 
W+G-sDe= 0 (16) 


mere G = kinetic energy input from the wind, 


D = dissipation within the mixed layer, and 
W = oie CS dz 
O 


mien substituting (16) into (15), expanding the total deriva- 


Meves, and solving for oT /ot yields 


Ze 





dT _ 2Q4 2(G-D dul. dvT. owl 


it Seema = ae a“ 


O'W Pp 0 
which is the predictive equation for change in the mixed- 
Mayer temperature due to convection, entrainment mixing, or 
advection. Setting Q, = 0, and/or (G-D) = 0, and/or Ties; 
and w = 0 allows isolation of the cooling mechanisms for use 
mecividually and in combination. Substitution of (17) into 
(14) yields 


eam Ra + _2(G-D) : vie ; ae +w (18) 
ReePueils ~ yp BCT, - TD aX dy 





fred 15 the predictive equation for changes in the mixed- 
layer depth as used in the present model. 

A partitioning procedure was used to divide the total 
available stress, ae from the boundary layer into one frac- 


meen, tor the production of turbulence, To? and the remainder 


cc. - ees hemmcicmecneration Of wind-driven currents. 
h 
Um 2 hn dee Cz ) (19) 
where ‘a Potmestwcesseavartlable from boundary-layer model, 
hy aeitee dt aver dept, 
zy = roughness length, and 
te-eameonstant (10) (after Fraim 1973). 


Kraus and Turner (1967) assuming a value for TO: 
Specified the mechanical energy input (G) in equations (16), 


7), and (18) as 





G = on) aie (20) 





Sravatational acceleration, 


= 

=f 
O 
a 
‘a 
QQ 
I 


a = (1/p,) (do, /dt), coefficient of expansion 


on air density, 

Pp, = sea water Gems ity, 

U, = atmospheric friction velocity, = ee: and 

fe weadt part Of atmospheric surface stress used to 


PROGguCce “EUAbDULenNt mi xine . 
The dissipation (D) was set equal to zero. 

The calculation for the ocean velocities was by means 
of Ekman-spiral theory. The theory specifies a depth (Z,) 
meewnich the Current direction has reversed from its surface 
@arection, the Current is EXP(-7) of its surface value, and 
Below which the wind has no direct influence. The velocity 
@omponents of the ocean within the spiral layer are 


(tT, -T,)cosa 
u = ———  EXP(-BZ)[cos(-BZ)-sin(-BZ)] (21) 


Hy % 
ee an? 
W 


(T,-T,) sina 
v= —+———_ EXP(-BZ) [cos(-BZ)-sin(-BZ)] (22) 
3, (Cae) = 
W 


Where Z 1/B, specifies the Ekman depth, 

B= [£/2u/o,1%, 

f = coriolis parameter, 
Oo = sea water density, 

Dene Cmte dimcddy coetficient of viscosity, 


z = depth, and 


@ = atmospheric inflow angle. 


Di 





Integration of (21) and (22) over the Ekman depth (2.5) gives 


the net mass transport in the Ekman layer 


Me ee - T Jcosa/t (2:5)) 


M 
y 


Cis ~ tT Jsina/f. (24) 
From these mass transport equations it 1s possible to calculate 


the mean ocean velocities in the mixed layer 


uZ, = M,/o,, (25) 


vz, = Mf P,,: (26) 


The horizontal velocities in the layer are related to the ver- 


tical velocity w by the continuity equation 
OW du dV 
see - (SE + SS). (27) 


Equations (17) and (18) give two equations for three 
unknowns (T.,h,,T))- A return flow is calculated for the 
layer between 130 meters and the new mixed-layer depth. It 
is assumed to have the same mass transport as the opposing 
Current in the mixed layer, and to be centered on a point 
2/3 of the distance from 130 meters to h,: The below-layer 
gradient is adjusted by the return flow, upwelling, and down- 
welling. The new i is the temperature at the intersection 
of the new below-layer gradient and the new mixed-layer 


Mm pth. this relation is given by 
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where dT/dz = is the below-layer gradient, 


Tp = temperature at top of the undisturbed layer 
= 20 C, and 
Zp» = depth of top of undistirbued layer = 130 m. 
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Bei locUSSION OF RESULTS 


It was the intent of this study to use a sea-air inter- 
acting numerical model to simulate the effects upon the 
@eean OL a hurricane moving across it. The reasons for con- 
aMie@eine the study were: 1) to investigate the cause of low 
sea-surface temperatures in the wake of a hurricane, and 2) 
memcetermine the effects on the ocean for different speeds 
of storm movement. 

Changes to the mixed-layer depth are controlled by Equa- 
Mrons (17) and (18). Convective mixing 1s represented by 
the first term in both equations. Entrainment mixing in- 
moves the second term in both equations. Advection is 
described by the remaining terms in the equations. Accord- 
moto the two equations entrainment mixing and convection 
can only cool and deepen the mixed layer. It has been as- 
sumed that Qa 1S positive by assuming a wind-dominated 
Mmeorne. On the other hand, since ocean velocity values can 
Memo megative and positive, the advection terms can both 
warm and cool, shallow and deepen the mixed layer. 

Hiomaiscussion 1s divided into three parts: 1) a com- 
parison with the results of Fraim (1973), 2) a discussion of 
the relative PPO Geavcesimetniis model] of the mechanisms dis- 
cussed above, and 3) a discussion of the effect of the hurri- 


Mane translational speed on changes in the ocean. 
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A. COMPARISON WITH FRAIM'S RESULTS 

Before attempting to move the hurricane, it was necessary 
to establish that the stationary version of the present model 
yielded reasonable results. For the initial experiment, the 
moo x 5 grid (Figure 1) was defined with a 3 km grid spacing 
as was the 100 x 1 grid in Fraim's model. Profiles of values 
@t JZ hours for the present model are shown in Figure 2. A 
comparison of maximum and minimum values for Fraim's model 
and the present model are listed in Table I. 

The first point of comparison between the two studies is 
Mone hurricane portion of the model. The initial conditions 
for both hurricanes were the same: air temperature i = 29C, 
radius of maximum wind ee 27km, and maximum tangential wind 
Ve. = 30m/sec. There were some differences. Fraim did not 
ee train the maximum tangential wind or the radius of maxi- 
mam wind. In the present model, for this experiment, the 
mMaXimum tangential wind was forced to increase sinsuoidally 
from 30m/sec to S50m/sec at 12 hours. Also, a limit of 
i < 60 km was set on the radius of maximum wind. 

iMiemeancential Wind protile (bottom, Figure 2) 1s char- 
acterized by a rapid increase with radius to its maximum 
value i mii EQmatrone2)), a slow decline to its 
value at Ve (Equation (1)), and a more rapid decline to the 
outer boundary value. From i to 2s the wind is allowed to 
mop off to Sm/sec. At ie the ¢hange in the wind profile is 
abrupt and therefore easily detected. The radius Ls tends 
meoe the limit of applicability of the asoumptclon of axi- 


Symmetry at upper levels, although the assumption is good to 
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faible i.) COlparison Between Some 12-Hour Values of the 
Fraim Radial Stationary Hurricane Model and the 
Mrametrical Ootataionary Hurricane Model. 


Fraim Model Diametrical Model 

Max ORV 48 as 
Max DW 05 nova 
Max UW S55 1.88 
Min MLD 4.31 Zale 2 5 
Max MLD Sli 55.20 
Min MLT (IAS 28.84 
Max MLT 30.43 30.00 
Max TW BS ebite 50 
RMW 42 60 
Area average ocean 

heat loss 930 cya 
Many, TT 19.80 26.90 
Max TT 29.99 30 
ORV - ocean radial velocity (m/sec) 


DW, UW - downwelling, upwelling (m/hr) 


MLD lec Gaver deme (im) 

MLT Ditccd tivewotemperature (C) 
TW - tangential wind (m/sec) 

RMW - radius of maximum wind (km) 
iy - thermocline temperature (C) 


Heat loss (cal/cm’) 
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Georeater radius at lewer levels. Because i expands with 
intensification (Tr, = 225 km at 12 hours vs. 150 km at 0 
Moeurs), the wind velocity value at the outer boundary of 
the hurricane model has increased from S5m/sec to 12.5m/sec. 
The maximum tangential wind (Table I) was about 10% greater 
imerraim's model than that of the present model due to lack 
of constraint on the maximum wind value. 

The ocean radial velocity (actually the x-component of 
the horizontal velocity) (middle, Figure 2) responds directly 
to the stress imparted by the tangential wind (Equations 
(21) and (22)). Ocean velocity values are calculated by 
imieecrating the x-component of the Ekman current over the 
maxed-layer depth. Fraim integrated the mass transport over 
meonmallower depth thus giving a larger velocity (Fraim's 
mesults, Table I[). 

As the ocean radial velocity directly follows the tan- 
foetal wind, the ocean vertical velocity (top, Figure 2) 
jmellews the gradient of the radial velocity. The result is 
that the maximum upwelling occurs at the point of maximum 
mnease wot Ocean radial velocity, 1.e€., between the center 
of the storm and the radius of.maximum wind. Maximum down- 
welling occurs at Le where the ocean radial velocity drops 
Sharply. Fraim's model upwelling (Table I) is stronger be- 
cause the radius of maximum wind was smaller. The Fraim 
model had no major downwelling because tangential wind pro- 
file was not altered in the r, to 2X0 region. 

The radius of the maximum mixed layer depth, minimum 


mixed-layer temperature, and minimum thermocline temperature, 
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(Figure 3) is the same, and coincides with the radius of 
maximum wind. Coinciding with the radius of maximum upwell- 
ing are the minimum mixed layer depth, and pips of warming 
im the mixed-layer ) and thermocline (T)) Polperatures. 
ime pips appear to be due to a lack of cooling due to upwell- 
ing. The minimum mixed-layer temperature (Table I) is much 
smaller in Fraim's model because it did not have the thermo- 
ominme Layer return flow discussed in Chapter II B. 1. Lack 
of a return flow of water allowed strong cooling of the 
thermocline layer in the Fraim model which created a large 
T,-T, difference (Table I) and retarded deepening of the 
feexeca Jayer. Consequently, the area-averaged heat loss in 
Poaim's model was more than twice the loss in the present 
model (see Table I}. Most of the heat loss is concentrated 
tomemam 100 km of the center. Advection of water both hori- 
womeally and vertically acts to cool or warm the thermo- 
memme layer. Warming in the upper portion of this layer 
increases the Doerr yeonecraarent while cooling decreases 
miemoradient. A change in the below-layer gradient is thus 
meagelatea with the temperature at the top of the thermo- 
mmc, that 1S, the thermocline temperature, since the tem- 
Merature at the level of no-motion is held constant. 
Stronger upwelling in Fraim's model (Table I) accounts for 
mae shallower mixed-layer depth minimum and seems to indicate a 
higher maximum mixed-layer temperature. The comparison be- 
tween the models shows the results are similar, however, due 
to the differences mentioned above Ss present model yielded 


more. reasonable results. 
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B. DISCUSSION OF COOLING MECHANISMS BENEATH A MOVING STORM 

As the comparison of the present model with Fraim's 
Was considered satisfactory, it was decided that the results 
warranted progressing to a moving storm. The storm was al- 
iewead tO interact with the ocean until a maximum wind of 
50 m/sec was achieved so as to allow oscamemnaveon of the 
Seca Seresponse to the Same storm. A new limit of eS 
km was placed on the radius of maximum wind. The 200 x 3 
grid spacing was expanded to 6 km giving a 12 km by 1200 km 
Mameanole in themopen oG€ean. Starting point for the center 
of the storm was the 900 km point (Figure 8). Movement of 
eiewstOnrm Was from right to left at 1 grid increment per 2, 
moe ana / time @ncrements, or 6 km per 16, 24, 40, and 56 
fomuMces Wikch corresponds to 11.84, 7.89, 4.75, and 3.39 
knots respectively. 7 

PetunescrOs>-scetion Of the 4.75 knot storm was drawn 
moreethne mixed-layer temperature response (Figure 4, Table II) 
and mixed-layer depth response (Figure 5, Table III) of the 
meeanm to the mechanisms indicated in the figures. The 750 
Km point was chosen for this section because during the 30 
Nour forecast period it experienced passage of the most in- 
Memestilp parts of the storm i.e. the radii of maximum wind, 
mtesradii Of maximum upwelling, and the no-stress zone of 
macy e. Figure € can be used to see the relative locations 
mon these none: 

As shown in Figure 4; advection alone has virtually no 
impact on the mixed-layer temperature, though in the absence 


of deepening (Figure 5) shows a slight warming as mentioned 
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my ChapewtIi. A. Entrainment mpxing and convection appear 
to dominate outside the radii of maximum upwelling with the 
maximum effect being inside the radii of maximum wind. Of 
the two, entrainment mixing is the stronger cooling mech- 
emism. It is interesting to note that convection alone con- 
tributes to the small amount of cooling within the radii of 
maximum upwelling. All mechanisms combined appear to con- 
tribute to cooling between the radii of maximum upwelling, 
although no changes occur during passage of the eye. It ap- 
pears that horizontal ocean velocities have little direct 
effect on the mixed layer, but that vertical velocities have 
Peoenificant effect. 

As shown in Figure 5, advection alone has a small effect 
mimenangcing the mixed layer. Again, entrainment mixing and 
convection appear to dominate outside the radii of maximum 
Mewelling with the maximum effect being inside the radii of 
maximum wind. Of the two, entrainment is definitely dominant 
mn deepening the mixed layer, although inside the radii of 
maximum upwelling convection is the only mechanism that 
deepens the mixed layer. Inside the radii of maximum upwell- 
ing, advection alone is responsible for shallowing the mixed 
Bayer. In conjunction with entrainment mixing and convection, 
the effect of advection increases. 

One cannot completely separate the mechanisms since there 
1s a complex interaction between them in the ocean. However, 
Leipper (1967) interpreted baththermograph data before and 
after passage of the hurricane Hilda in the Gulf of Mexico 


during 1964 as follows: 1) cold water was upwelled due to 
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divergence in the region near the center of the storm, and 
2) warm water farthest from the storm path gave evidence of 
having been mixed by surface cooling or by mechanical mixing 
or both. As shown a Fieure 9yetmne lowest mixed-layer tem- 
peratures are in the vicinity of the storm Senter The 
warmer water toward the boundary of the storm (Figure 8) 
shows evidence of mixing (implied by mixed layer deepening). 

A set of profiles of temperature change from the initial 
temperatures in the ocean due to hurricane passage is shown 
in Figure 6. The profiles show cooling down to 30 m (except 
the profile for advection only) then warming down to the 
mixed-layer depth. Entrainment of colder thermocline water 
cools the mixed layer and warms the thermocline layer to 
mormea mew mixed layer. At the mixed-layer depth in the 
model there is a discontinuity denoted earlier by le - T)- 
Below the mixed layer only those mechanisms involving advec- 
mone snow distinct cooling. Return flow in the thermocline 
layer and upwelling are responsible for the cooling. The 
750 km point is 120 km behind the storm center at 30 hours. 
AS a2 result Figure 6 agrees with Federov's (1973) Type A 
Miermal structure change profile of two layers of cooling 
separated by a layer of warming. 

The effect of the return flow in the thermocline layer 
mimealisSo be seen in Figure 7. Entrainment mixing and con- 
~vection Singly and together have an insignificant effect on 
the below-layer gradient. The rapid decrease of the below- 
mayer gradient for advection alone shows that for intense 


upwelling there is strong cooling in the upper part of the 
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thermocline layer indicating a more intense upwelling than 


for advection alone. 


C. EFFECT OF STORM MOVEMENT ON OCEAN RESPONSE 

The second purpose listed for the study was to determine 
the effect of movement of the storm on the ocean's response 
to the storm. In Table IV is listed some maximum and minimum 
values of key parameters for several storm translational 
speeds. 

Changes in the mixed layer due to storm movement are 
comparatively easy to explain. Slower-moving storms have 
femme Cime at a given location for entrainment mixing and con- 
fomeelon €O COOl and deepen the mixed layer. For the sta- 
tionary storm the downwelling maximum nearest the storm 
Smee, (see Figure 10) has added to the deepening but per- 
haps retarded the cooling. The anomalous spikes of down- 
welling (for instance at the 900 km point in Figure 10) are 
due to the SampeGtEoOpmiiethe tangential wind at the storm 
boundary and the mixed-layer depth spikes generated due to 
the stationarity of the storm. 

metreasing heat flux to the storm with increasing storm 
speed (Table IV) is due to fast moving storms extracting more 
meau because of continuously encountering uncooled 30 C 
water. Radial heat flux (Table IV) is due to heat flow in 
Mimrents OUt through the mixed layer and in through the 
thermocline layer. For a stationary storm advection supplies 
more heat than it removes. This seems to be the source of 


Warm water for the mixed-layer temperature pips (Figure 11). 
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Table IV. 


max ORV 

max DW 

max UW 

min MLD 

max MLD 

min MLT 

max MLT 

area average 
ocean heat 


loss 


radial heat 
aX 


meat lost to 
storm 


*&see Table I for abbreviations 
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The apparent reversal of radial heat flux values for the two 
fastest storms in unexplained. 

A sustained upwelling for 18 hours duration would account 
for the shallow minimum mixed-layer depth (Table IV). By the 
Pauation for continuity (Equation (27)) the integrated ver- 
tical velocity is equal to the horizontal velocity divergence 
memes the integrated depth. Therefore the maximum upwelling 
value for the stationary storm is small because of the mixed- 
layer depth in which it occurs (Figures 10 and 11). Fora 
moving storm the leading upwelling zone the upwelling is more 
intense because it occurs in the deepest mixed layer (Figures 
8 and 9). Maximum downwelling intensities are governed by 
meeeed-layer depth also (Figures 8, 9, 10 and 11). Maximum 
Values for the ocean radial velocity are much the same be- 
Cause the mean mixed-layer depths at the radius of maximum 
ocean radial velocity are nearly the same. The ocean radial 
velocity profiles are altered slightly due to mixed-layer 


depth considerations (Figures 8 and 10). 
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IV. SUMMARY AND CONCLUSION 


It was the intent of this study to use a numerical model 
Porsimulate the effects of a hurricane moving across the 
feean. In such numerical models the inherent limitations in- 
clude the finite differencing scheme as well as the initial 
and boundary conditions. Imposed limitations involve assump- 
tions associated with applying the model. 

For this study a simple model was used for both the hur- 
ricane and the ocean. A bound was set on the hurricane 
maximum wind in order to standardize the storm input to the 
ocean. There were several key assumptions about the way in 
which the ocean and storm interact: 1) wind stress was par- 
Mmieeroned into making currents and turbulence based on the 
mexecd-layer depth, 2) ocean currents were assumed to be in 
Ekman balance and to respond immediately to changes in the 
Meness Pattern, 3) it was assumed that there were no gradi- 
pes) OL temperature or ocean velocity across the path of the 
Storm, 4) the storm was assumed to be wind-dominated, and 
5) since the storm was mature at initiation the total stress 
imparted exponentially until 95% of the total stress was 
imparted to the ocean at 9 hours. 

In answer to the two main questions of interest: 1) 
meomine Of the ocean appears to be primarily due to entrain- 
Ment and convection except in intense upwelling within the 
heeton OL Maximum wind. In this region all mechanisms con- 


merpute to cooling of the mixed layer, 2) slow moving storms 
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cool and deepen the mixed layer more than do fast moving 
storms. Waves in the mixed-layer depth due to internal 
waves on the mixed-layer thermocline layer interface postu- 
mcd py Gerster (1970) were not seen. This is due to the 
thermodynamic nature of the model under study here compared 
to the hydrodynamic nature of Geisler's model. 

In summary, the results are satisfactory within the limits 
of the model. The logical course to follow is to reduce the 
number of assumptions used. In a model with as many simul- 
mmeous ifteractions as this one, the task is one to be taken 


slowly, in small increments. 
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